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We study the global analysis for parton distributions as a function of the QCD strong 
coupling strength and present a new series of distributions that span the range 0.110 < 
O-sijnz) < 0.128. We use these distributions to explore the correlation between and 
the gluon distribution; the viability of global analysis as a method to measure as] and the 
dependence on ag of predictions for W, Z, inclusive jet, and Higgs boson production {bb —> H 
and gg —> H) cross sections at the Tevatron and the LHC. We hnd that the uncertainty in 
as is the dominant source of uncertainty for inclusive jet production at moderately small pt 
and for Higgs production at intermediate masses in the standard model. 
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1 Introduction 

The theory of Quantum Chromodynamics (QCD) depends on the fundamental gauge cou¬ 
pling strength and suitably defined quark mass parameters ruf. For applications to hard 
scattering processes with hadrons in the initial state, we also need the universal parton dis¬ 
tribution functions (PDFs) that characterize the partonic structure of those hadrons. These 
PDFs are determined by global QCD analysis, using input from a variety of well-established 
experimental measurements. 

In the standard CTEQ analysis [1], the focus is on determining the parton distribution 
functions. The value of as{mz) is fixed at the world average value, which is dominantly based 
on dedicated measurements such as those from LEP, where it was determined in processes 
that are free of the complications of hadronic structure. 

However, the interplay between the coupling strength and strong dynamics is an 
interesting subject in itself. Many attempts have been made to extract as{mz) from indi¬ 
vidual experiments at the HERA ep collider, at the Tevatron p] collider, and from combined 
analyses of several hadronic experiments. For such studies, as well as to assess the additional 
uncertainty in predictions caused by the uncertainty in as{mz), it is important to have PDF 
sets available that are based on a range of different values for ag- The purpose of this paper 
is to fill that need. 

We provide here a series of PDF sets that span a range of coupling strengths from 
O-sijnz) = 0.110 to 0.128. These PDFs extend and update the CTEQ6 global analysis 
[1]. They are successors to the series that was determined in the CTEQ4 analysis [2]. 
Complementary to the physics probed in [1], we will discuss several issues related to the 
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dependence on a*: the correlation between Og and the parton distributions—particularly 
the gluon distribution; the viability of using global analysis of hadronic processes to measure 
as] and the dependence on of physics predictions for W, Z, jet, and Higgs boson cross 
sections at the Tevatron and the Large Hadron Collider LHC. 

We provide two different sets of hts—CTEQ 6 A and CTEQ 6 B—corresponding to two 
different common dehnitions for as{fi) that are dehned in the next Section. 

2 Definitions for 

The dependence of the QCD coupling strength as{fi) on the momentum scale /r is gov¬ 
erned by the renormalization group equation (RGE), which in next-to-leading order (NLO) 
perturbation theory is 

jjida/djji = cio^ -\- c^c? , (1) 

where Ci = —I3q/2t[ with = 11 — (2/3)nj, and C 2 = —/ 9 i/ 87 r^ with fdi = 102 — (38/3)nj. 
The coefficients are functions of the number of active quark ffavors Uf. These coefficients 
change discontinuously as n passes the mass m/ of each quark ffavor and the integer n/ 
jumps up by 1; nevertheless at NLO, the function as{fi) is a continuous function of /i at 
these thresholds. The solution to Eq. (1) therefore depends on a single integration constant, 
which is the fundamental coupling strength parameter of QCD that is usually chosen to be 
as{mz). 

Since Eq. (1) is truncated at 0{a^), there are inhnitely many dehnitions of as{fi) that are 
formally equivalent because they differ only in higher orders of the perturbative expansion. 
Two choices that are often used in QCD phenomenology are: 

• Def. A. The original NLO dehnition [3] is given by 

as(h) = C3 [1 - C4ln(L)/L]/L , (2) 

where L = ln(yU^/A^), C 3 = —2/ci, and C 4 = —2 02 / 0 ^^. The parameter A depends on 
nj, and hence takes on different values A„^ when /r crosses each quark mass threshold. 
Previous CTEQ global analyses have all used this dehnition for the running coupling. 

• Def. B. An alternative is to solve the truncated RGE (1) exactly. The publicly avail¬ 
able evolution program QCDNUM [4], and many HERA analyses use this dehnition. 

We have shown previously [5] that these two forms are numerically quite similar in the 
region Q > 2 GeV where we ht data. In the following, we present results based on both 
dehnitions,^ 

^The form used for PDF fitting by MRST [6] is different from both of these definitions, but is numerically 
very close to Def. B [5]. The definition used by the Particle Data Group [7] lies between Def. A and Def. B. 
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3 The CTEQ6 series of global fits 


To study the interplay between the strong coupling strength Og and the PDFs, one can 
either perform a fully global QCD £t by varying as and the PDF parameters simultaneously 
to examine the neighborhood of the overall global minimum in or one can perform a 
series of fits to the PDF parameters at various hxed values for We have explored both of 
these approaches, but concentrate here on the second approach since it is most convenient for 
general collider physics applications. We therefore present here a series of PDFs—the CTEQ6 
Q!s-series. Some of their physical implications are discussed in Sec. 6. Earlier PDF hts with 
a range of as{mz) were obtained in [2,8]. Fits with a range of as{mz) and their interplay 
with the gluon distribution were studied in connection with the possibility of coevolution 
with light gluinos in [9]. 

The PDFs in the new series were obtained for 10 values of as{mz) = 0.110, 0.112, ..., 
0.128, where mz = 91.188 GeV is assumed. The theoretical assumptions and functional 
parametrization of these PDFs at the initial momentum scale /io = 1-3 GeV are the same 
as in the previous GTEQ 6 analysis [1,10]. The experimental input is slightly updated.^ 
These PDF sets are designated as GTEQ6A110, ..., GTEQ6A128 using Def. A for and 
GTEQ6B110, ..., GTEQ6B128 using Def. B for as- The hts GTEQ6A118 and GTEQ6B118, 
which have as{mz) = 0.118, are nearly the same as GTEQ 6 M or GTEQ6.1M, but are not 
identical to either due to minor updates in experimental input, and in the case of GTEQ6B118 
the different dehnition for as- All 20 of these PDF sets will give very similar physical 
predictions for most applications. 

Figure 1 shows the quality of the global hts, as measured by the overall fhe ht to 

~2000 data points, as a function of as{mz)-^ The two curves, both approximately parabolic, 
are smooth interpolations of the above series of hts. The curves are very similar, with the 
Def. A curve being slightly narrower and slightly farther to the left because Def. A has a 
little more rapid variation of as{fi) with /r. The minima of these curves are at as{mz) = 
0.1172 and 0.1176, close to the current world average (0.1187 ± 0.0020) [7] and very close to 
the value 0.1180 used in GTEQ 6 M and GTEQ6.1M. This similarity in as is an impressive 
demonstration of consistency between QGD theory and experiment, since the global QGD 
analysis is based on hard scattering data with hadronic initial states, while the determination 
of the world average value for comes mainly from totally diherent physical processes such 
as annihilation, r-decay, and even lattice gauge theory calculations with quarkonium 
spectra as input. 

Note that the range of as{mz) covered by the GTEQ 6 AB series is much wider than 

^The input experimental data consist of BCDMS (muon F 2 on hydrogen and deuterium) [11], HI (e^ 
F 2 ) [12], ZEUS (e^ F 2 ) [13], NMC (muon F 2 on hydrogen and deuterium) [14], CCER (neutrino F 2 and 
xFs on iron) [15], E866 (Drell-Yan muon pairs from pp and pd) [16], CDF (W-lepton asymmetry) [17], CDF 
(inclusive jet) [18] and DO (inclusive jet) [19]. The update consists mainly of including the data from the 
first of the three references in [12] which was inadvertently omitted in CTEQ6. 

^This figure is an improved version of one appearing in [5]. 
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the currently accepted la error range of the world average quoted above. We explore this 
extended range because the lowest and highest values (0.110, 0.128) represent outlying val¬ 
ues that have been obtained by some individual experiments that are included in the world 
average. The hts with extreme values of asijnz) will be useful for some specialized applica¬ 
tions. In plots shown later in this paper, we show almost all of this range, but reduce it to 
0.110 to 0.126 to be symmetric about the CTEQ6 value of 0.118. Note that increases by 
~ 100 above its minimum at the extremes of this reduced range, which makes it consistent 
with a 90% conhdence range for the global £t according to results of our previous analyses. 
However, the reader must keep in mind that this range of as{mz) is larger by a factor of ~4 
than the uncertainty range corresponding to a “la” error band based on the world average 
data. 


4 The gluon distribution and 

The gluon distribution is strongly correlated with as in the global QCD analysis. This can 
be seen in Fig. 2(a), which shows the gluon distributions g{x,fi) from the CTs-series PDFs as 
a function of momentum fraction x at scale /i = 3.162 GeV {g? = lOGeV^). For clarity of 
display, the horizontal axis is scaled as x^^^ and the vertical axis is weighted by 

For comparison. Fig. 2(b) shows the uncertainty band (shaded area) of the gluon dis¬ 
tribution due to sources other than ag- This uncertainty was computed at a hxed value 
oisijnz) = 0.118 by the Hessian method [20], using the 40 eigenvector basis sets of GTEQ6.1 
[10]. The GTEQ6.1M (solid) and the new GTEQ6A118 (dashed) distributions, which are 
very similar, are also shown in Fig. 2(b). The two hgures are combined in Fig. 2(c), where, 
in order to highlight the differences, the results are shown as ratios to the GTEQ6.1M dis¬ 
tribution. 

There is a clear systematic trend in the cts-series for the gluon distribution function. 
Fits with larger as{mz) have a gluon component that is weaker at small x and stronger at 
large x. The behavior at small x results from the fact that every occurrence of g{x,fi) in a 
cross section formula is accompanied by a factor of as, so when as is made larger, the gluon 
distribution becomes smaller in order to maintain agreement with the large amount of data 
at small x. The behavior at large x, where the direct experimental constraints on the gluon 
are weaker, is dictated by the momentum sum rule: the total momentum fraction carried 
by gluon -|- quarks must be equal to 1, and the momentum carried by quarks is tightly 
constrained by DIS data. 

Another feature seen in Fig. 2 is that the gluon distributions for different as{mz) values 
all nearly intersect at a common value x ~ 0.1. This occurs simply because the func¬ 
tion g{x, n, asijnz)) varies rather slowly and smoothly with asimz)), so its dependence on 
(y-sifnz)) can be approximated rather well by the hrst order (linear) term of a Taylor series. 
This linearity can also be seen by the nearly equal spacing of the curves in Fig. 2(a) and (c). 

Gomparing the range of variation of gix, jj) due to the variation of as, to the uncertainty 
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range due to other sources of error, we see that the two are comparable throughout most of 
the domain in x. The combined uncertainty on the gluon distribution is therefore somewhat 
larger than the previously published uncertainties, which were obtained at hxed as{mz)- 

Figure 3(a) shows the gluon uncertainty at p = 100 GeV. At this larger momentum 
scale, the overall uncertainty is much smaller, and the as contribution to the uncertainty is 
generally smaller than the uncertainty due to the causes that are included in the Hessian 
analysis. Figure 3(b) similarly shows the uncertainty for the up quark distribution at /i = 
100 GeV. The uncertainty for the quark is much smaller than for gluon (note the different 
y-axis in the graph) and again the as contribution is small compared to the other sources of 
uncertainty. 

5 Can we determine cts from the global analysis? 

It is natural to try to determine the coupling strength from the global QGD ht to hadronic 
processes. It appears straightforward to do so, since one can simply treat the parameter 
(y.s{f^z) as one of the htting parameters. We see from the minima in Fig. 1 that the resulting 
“best £t” value of as{mz) is around 0.1174, which is very close to the world average value. 
The difficult question, however, is what uncertainty should be assigned to this measurement. 
The answer to that question determines whether this method is competitive with measure¬ 
ments that are independent of the complications of hadron partonic structure. The fact that 
the value of as{mz) in global analysis is strongly correlated with the rather uncertain gluon 
distribution, as discussed in the previous section, suggests that caution is needed. 

Referring again to Fig. 1, the range of as{mz) (horizontal axis) allowed by global analysis 
based on minimization of is set by the increase of (vertical axis) above the global 
minimum that we allow. We refer to the allowed increase of the tolerance, Ay^. 

Recent studies of uncertainties of PDFs by various global analysis groups [1,6,21] have 
concluded that a reasonable tolerance Ay^ must be rather large, in the range 50 - 100 for 
the ~2000 points in present-day data sets, to dehne an approximate 90% conhdence range. 
Making the specific choice Ay^ = 100, the corresponding range of as{mz) is from 0.1093 to 
0.1247. Assuming that range to be the 90% conhdence range for a gaussian distribution, it 
corresponds to as{mz) = 0.1170 ± 0.0047 for a “la” error range. Thus the measurement 
of as{mz) provided by PDF htting agrees very well with the Particle Data Group world 
average of 0.1187 ± 0.0020 [7], but has more than twice its uncertainty. Hence the PDF 
result cannot be used to substantially reduce the uncertainty in agimz) at the present time. 

We can gain some insight on how as is constrained in the global QGD analysis by 
examining the dependence on as{mz) of the y^ values for each individual experiment that 
contributes to that analysis. These y^ “parabolas” are shown in Fig. 4. The curves represent 
smooth interpolations of the results from the 10 hts in the GTEQA a^-series. The vertical 
axis in each graph is the y^ value per data point in the ht to that experiment, while the 
horizontal axis is as{mz)- 
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It is apparent in Fig. 4 that the sensitivities of the various experiments to Og vary 
greatly. Some of the curves are approximately parabolic with a minimum within the range 
probed, while others merely constrain the value of as{mz) from above or below. The global 
minimum seen in Fig. 1 is due to the combined constraints of all the experiments. Since 
different experiments prefer different values of as{mz), which are not always consistent with 
each other if strict statistical criteria (“Ay^ = 1”) are applied to each experiment, the global 
minimum represents a compromise that is difficult to interpret as a “measurement” in the 
traditional sense.^ In particular, there is no clear way to assign a statistically meaningful 
error to the measurement. Rather, the error is dominated by systematic effects that can 
only be estimated. 

6 Dependence of predicted cross sections on Og 

Here we present predictions for several important processes at the Tevatron (pp at a/s = 
1.96TeV) and the LHC (pp at a/s = 14TeV). The results show how the uncertainty in 
propagates to uncertainties in physical predictions. At the same time, the results show to 
what extent accurate measurements of these cross sections could be used to constrain 

We show predictions using the Def. A form for ag- Results for the Def. B form are 
very similar. The hgures show the variation of the predictions for the range of ag{mz) from 
0.110 to 0.126. We again remind the reader that this range is much larger than the actual 
uncertainty of as{mz)', i.e., the full variation of the physical predictions shown in each hgure 
extends beyond the actual uncertainty from ag. The shaded region in each hgure shows the 
range of uncertainty due to sources other than Og, as calculated from the eigenvector basis 
sets of CTEQ6.1 using the Hessian method [10,20]. These are assumed to estimate the 90% 
conhdence range. 

6.1 W and Z production 

Figure 5 shows the cross section for W~ production at the Tevatron.® production is 
identical except for y —>■ —y. The left plot shows da/dyw versus the W rapidity yw- The 
right plot shows the difference from the prediction of CTEQ6.1, which has ag{mz) = 0.118. 
The curves are the CTEQ6A cts-series, and the shaded band is the range of uncertainty for 
hxed ag{mz), calculated using the Hessian method. For this process, the Hessian uncertainty 
range is about ±5% of the central prediction. The variation with ag is smaller, on the order 
of ±2%, even for the extreme range of ag{mz) from 0.110 to 0.126. 

"‘We should point out that the curves shown in Fig. 4 are not to be compared directly to those obtained 
by individual experiments in their respective determinations of ag. The points on our curves correspond 
to values evaluated using constrained fits to the full global data set, not just to the data of a single 
experiment. 

®A11 cross sections shown for and production are given in nanobarns, with the leptonic decay 
branching fraction included. 
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Figure 6 shows the cross section da/dyz for production at the Tevatron. Again the 
Hessian uncertainty range is ~ ±5% and again the variation with as{mz) is ~ ±2% for the 
extreme range of 

Figure 7 shows the cross section for W~ production at the LHC. The Hessian uncertainty 
range is again of order ±5%, but the variation with as is larger than for the Tevatron, of 
order ±5% for the large range of asijnz) that is shown. Figure 8 shows the process of 
production at the LHC. It has a larger cross section at large rapidity (because u{x) > d{x) 
for the valence quarks) and a similar range of uncertainty in the prediction. The difference 
between the central dashed curve and a horizontal line in Figs. 7 and 8 shows the effect of 
updates in the htting between CTEQ6.1 and CTEQ6A118: the change is well within the 
estimated PDF errors. 

The cross sections for production at the LHC are closely tied to the gluon distri¬ 
bution, since the leading-order process in proton-proton collisions is u + d ^ W~^, which 
involves a sea quark; and sea quarks and the gluon are related at large y by the evolution 
equations. (Proton-antiproton collisions at the Tevatron are different because both u and d 
can be valence quarks—which causes the asymmetry in y that can be seen in Fig. 5.) Also, 
the next-to-leading order process u + g ^ d + involves an initial gluon directly. Hence 
the predictions for W production at the LHC are somewhat more sensitive to a* than are 
the predictions for the Tevatron. Nevertheless, the uncertainty associated with as at the 
LHC, for values of as that are consistent with the world average, remains small compared 
to the other PDF uncertainties as can be seen in Figs. 7(b) and 8(b). 

6.2 Inclusive Jets 

Figure 9 shows the a^-dependence of predictions for inclusive jet production at the Tevatron 
(in the CDF central rapidity region 0.1 < \y\ < 0.7), and at the LHC (in the region \y\ < 1). 
The as dependence is exhibited by plotting the ratio of the predicted cross section da/dpT 
to the prediction calculated using CTEQ6.1. As before, the shaded region is the uncertainty 
range for hxed as{mz) = 0.118 calculated from the eigenvector basis sets of the Hessian 
approach [10, 20] and the curves show asijnz) = 0.110, ..., 0.126. 

We observe that if we restrict as{mz) to a range ±0.003 which is the 90% conhdence 
range of the world average, the uncertainty due to as{mz) is small compared to the other 
PDF uncertainties at large pr- However, at moderately small pt the uncertainty in as{mz) 
adds considerably to the overall uncertainty for jet production. This comes about because 
the parton distributions do not depend very strongly on asimz) (see Fig. 3), but the hard 
cross sections do. 

6.3 Higgs boson cross section in SM and MSSM 

The uncertainty of the cross section predicted for Higgs boson production at the LHC as 
a function of Higgs mass is shown in Fig. 10. Figure 10(a) shows the uncertainty for the 
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gg —>■ H process, while Fig. 10(b) shows the uncertainty for the bb H mechanism. These 
cross sections were calculated at NLO using programs from [22,23] and [24] respectively. 
Both processes play an important role in Higgs physics: the gg process is dominant in the 
Standard Model (SM), while in Supersymmetry and some other generic extensions of SM, the 
bb ^ H process can be equally important or even dominant. The study of PDF uncertainties 
is therefore important for both mechanisms of Higgs production. 

The uncertainties based on the 40 eigenvector sets of CTEQ6.1 are shown as the shaded 
regions. These have been shown previously for gg —>■ H [25] and for bb H [26]. The 9 
curves in each figure show the predictions of the PDF sets with as{mz) = 0.110,..., 0.126 
in steps of 0.002 relative to the central value 0.118, which therefore corresponds to the 
horizontal line at ratio 1.0. 

Figure 10(a) shows that the uncertainty in substantially increases the uncertainty 
in the prediction for gg — H. Particularly for Mh — 300 GeV where the non-Os PDF 
uncertainty is a minimum, the uncertainty due to as is larger than the non-Os uncertainty. 
The strong sensitivity of gg ^ H to is of course not surprising in view of the af dependence 
coming from the leading order triangle diagram; indeed because the K-factor is large and 
positive for this process, the NLO corrections make the dependence on as even stronger.® 

Figure 10(b) shows that the contribution of the uncertainty in as to bb —>■ H production 
is comparable to the other PDF uncertainties for that process, so its effect is to add just 
a modest increase in the overall uncertainty for that process. This is not surprising, since 
bb ^ H has no direct dependence on ag at leading order, and the quark distributions at this 
momentum scale do not vary rapidly with as as seen in Fig. 10(b). 

Figure 11 summarizes the results of Fig. 10: The dotted curves are the PDF uncertainty 
at hxed (identical to the boundary of the shaded regions in Fig. 10) which are intended 
to show a 90% confidence range. The dashed curve is the uncertainty due to Os, calculated 
by interpolation for an uncertainty of ±1.64a in as{mz), where a = 0.002 from the world 
average and the factor 1.64 corresponds to a 90% conhdence range. The solid curve shows 
the combined uncertainty obtained by adding the two contributions in quadrature.^ The 
additional uncertainty due to a* is seen to be substantial for the gg —>■ H process. 

®This simple notion agrees quite well quantitatively with the results in Fig. 10(a). For instance, at 
Mh = 100 GeV, raising Us by 5% from 0.118 to 0.124 would be expected to raise the cross section by about 
13%—halfway between 10% (for a^) and 16% (for a^). But meanwhile, Mh = 100 GeV requires x « 0.007, 
where the increase in Us causes g{x) to decrease by about 3% according to Fig. 3(a), which would lower the 
cross section for gg ^ H hy 6%. Combining these two effects, the net change is an increase of 7%, which 
agrees well with the actual increase of 8% that is seen in Fig. 10(a). 

^Adding the errors in quadrature, i.e., treating the additional source of error due to as{mz) as independent 
of the other PDF errors, is the correct approach according to the Hessian approximation. For in that 
approximation, = Xo + where the Zi are linear combinations of the PDF shape parameters Ai'. 

Zi = J2j TijAj. When a new parameter such as {as{mz) — 0.118) is added to the set of fitting parameters, 
Xo becomes a quadratic function of that parameter, but there is no change in T^- because in the Hessian 
approximation one drops all contributions to x^ that are higher order than quadratic. 



7 Conclusion 


Previous CTEQ6 parton distributions were extracted from experiment assuming as{mz) = 
0.118, based on the world average value. The PDFs presented here were extracted using a 
range of alternative assumptions for as{mz), using a similar set of experiments with only 
minor updates. 

The new PDFs are named CTFQ6A110, ..., CTFQ6A128, and CTFQ6B110, ..., 
CTFQ6B128, where the “A” or “B” label indicates the choice of functional form for «<*(/!) 
(see Sec. 2) and “110” e.g. indicates as(91.188 GeV) = 0.110. Fortran programs to calculate 
these PDFs are available at the LHAPDF archive http://www-spires.dur.ac.uk/HFPDATA/. 

These PDFs can be used to hnd the range of uncertainty in predictions due to the 
uncertainty in as{mz). We hnd that the previous hts with as{mz) = 0.118 are adequate 
for most processes, because the uncertainty associated with as{mz) is smaller than the 
other sources of PDF uncertainty. However, the as{mz) uncertainty is important for some 
predictions that are particularly sensitive to a* or to the gluon distribution, such as inclusive 
jet production at relatively small pt (see Fig. 9) and Higgs boson production by the gg —>■ H 
process in the standard model (see Fig. 11(a)). This comes about because the intrinsic 
that is present in the hard scattering processes is only partly compensated by changes in the 
PDFs with as{mz) such as those shown in Fig. 3. 

These PDFs can also be used to study the constraints from global htting or to study 
individual experiments on the value of as{mz)- In the latter case, the new alpha-series 
PDFs allow one to take into account the strong correlation between as{mz) and the gluon 
distribution that is present at small momentum scales, as shown in Fig. 2. 
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Figure 1: The overall goodness-of-fit measure for global fits vs. as{mz), using Def. A 
(solid curve) and Def. B (dashed curve). 
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Figure 2: (a) Gluon distributions for CTEQ6A110 (short dash), CTEQ6A126 (long 
dash). These have asijnz) = 0.110 (short dash), ..., 0.126 (long dash) using Def. A for ag] 
(b) Gluon distributions with as{mz) = 0.118: GTEQ6A118 (dashed), GTEQ6.1M (solid), 
and uncertainty band from GTEQ6.1 eigenvector sets (shaded); (c) Both plots combined in 
ratio form. 
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Figure 3: Uncertainties at scale fi = 100 GeV for Gluon (a), and u-quark (b). 
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Figure 4: The values per data point for the individual experiments that are included in 
the global analysis, as a function of as{mz)- The number of data points is indicated on each 
hgure. 
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Figure 5: (a) Cross section da/dyw for W~ production at the Tevatron. (b) Same cross 
section with the CTEQ6.1 prediction subtracted. The curves are for as{mz) = 0.110 (short 
dash), ..., 0.126 (long dash) as in Fig. 2. 
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Figure 6: (a) Production of at the Tevatron; (b) Same process with the CTEQ6.1 pre¬ 
diction subtracted. 
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Figure 7: (a) Production of W at the LHC; (b) Same process with the CTEQ6.1 prediction 
subtracted. 
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Figure 8: (a) Production of IV~^ at the LHC. (b) Same process with the CTEQ6.1 prediction 
subtracted. 
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Figure 9: Uncertainties of inclusive jet predictions for Tevatron (a); and LHC (b). 
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Figure 10: Uncertainty of predictions for Higgs boson production at LHC for the SM and 
MSSM process gg ^ H via top triangle diagram (left); and the MSSM process bb ^ H 
(right). Curves (short dash to long dash) show the ratio (jff(as(m^))/(jff(0.118) for as{mz) = 
0.110,..., 0.126. Shaded region shows the uncertainty at as{mz) = 0.118 from other sources. 
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Figure 11: Percentage uncertainty of predictions for Higgs boson production at LHC at 90% 
confidence. Dashed curves are due to asijnz) uncertainty; Dotted curves are due to PDF 
uncertainty at fixed asijnz)] Solid curves are the combined uncertainty. 
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